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C. LANDIM, M. SUED AND G. VALLE 
Dedicated to Jozsef Fritz on his sixtieth birthday. 

Abstract. We consider the asymmetric exclusion process. We start from a 
profile which is constant along the drift direction and prove that the den- 
sity profile, under a diffusive rescaling of time, converges to the solution of a 
parabolic equation. 



1. Introduction 

Consider the asymmetric exclusion process evolving on the lattice "L d . This dy- 
namics can be informally described as follows : fix a translation invariant transition 
probability p(x, y) = p(0, y — x) = p(y — x). Each particle, independently from the 
others, waits a mean one exponential time, at the end of which being at x it chooses 
the site x + y with probability p(y). If the chosen site is unoccupied, the particle 
jumps, otherwise it stays where it is. In both cases, after its attempt, the particle 
waits a new mean one exponential time. 

The configurations of the state space {0, 1} Z are denoted by the Greek letter rj 
so that, for x in Z d , rj(x) is equal to 1 or 0, whether site x is occupied or not. For 
each density < a < 1, the Bernoulli product measure with parameter a, denoted 
by u a , is invariant. 

The macroscopic evolution of the process under Euler rescaling is described [H] 
by the first order quasilinear hyperbolic equation 

d tP + q-VF(p) = 0, (1.1) 

where F(a) — a(l — a) and q e M d is the mean drift of each particle : q = 
^2 z zp(z). Assume that the system starts from a product measure with slowly 
varying density po(eu). Under Euler scaling (times of order te _1 ) the density has 
still a slowly varying profile X e (t, su) which converges weakly (in fact pointwisely 
at every continuity point, 0) to the entropy solution of equation (|1.1|) with initial 
data pq. 

In the context of asymmetric interacting particle systems the Navier-Stokes 
equations takes the form 

d tP " + q-VF(p s ) = eY J du i (a i , i {p s )d UjP s ) , (1.2) 

where a is a diffusion coefficient. Three different interpretations have been proposed 
for the Navier-Stokes corrections : 

(a) The incompressible limit (|3], @]) : Consider a small perturbation of a 
constant profile ao : p% = ao + Assuming that this form persists at latter 
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times (p £ (t,u) = ao + eip(t,u)) we obtain from (|1.2J) the following equation for 

ife = (^(te _1 ,u) 

d m + e- 1 F'(a )q-cp E + (l/2)F"(a )q ■ = Oyfoo)]]^,^^ + 0(e). 

A Galilean transformation m £ (t,u) = ip s (t,u + e~ 1 tF'(ao)q) permits to remove 
the diverging term of the last differential equation and to get a limit equation for 
m = lim e ^o Trie 

d t m + (l/2)F"(a )g- Vm 2 = a id {a ) XX,« 3 m ■ 

(b) First order correction to the hydrodynamic equation ( 2 , 8 ) : Fix 

a smooth profile po : ^ ~~* an d consider a process starting from a product 
measure with slowly varying density po(eu). We have seen that under Euler scaling 
the density is still a slowly varying profile X £ (t, eu) which converges weakly to the 
entropy solution of equation Hl.lfl with initial data po- This second interpretation 
asserts that the solution of equation ljl.2|l with initial profile po approximates A e 
up to the order e : 

e-^-ff) - 

in a weak sense as s J. 0. 

(c) Long time behaviour ([5], [Q) : The third interpretation consists in analyzing 
the behaviour of the solution of equation l|1.2|l in time scales of order te~ l . Let 
b e (t,u) = p(te _1 ,u). From Q1.2f> we obtain the following equation for b e : 

d t b £ + £ -VVF(6 e ) - J2 d ^( ai ^ 9 ^) ■ 

To eliminate the diverging term e~ l q ■ VF(6 £ ), assume that the initial data (and 
therefore the solution at any fixed time) is constant along the drift direction : 
q ■ Vpo = 0. In this case we get the parabolic equation 

dtb e = ^d Ul (ai^(b e )d Uj b^j 

which describes the evolution of the system in the hyperplane orthogonal to the 
drift. 

Notice that while the first and the third interpretation concern the behaviour of 
the system under diffusive rescaling, the second one is a statement on the process 
under Euler rescaling. Interpretation (a) and (b) have been proved 0], |S] for 
asymmetric simple exclusion processes in dimensions d > 3 and a double variational 
formula for the diffusion coefficient was deduced. As one would expect, the diffusion 
coefficients of the two interpretations are the same and may be expressed by a 
Green-Kubo formula [T5|, It was also proved (Corollary 6.2, [S]) that the diffusion 
coefficient is strictly bounded below in the matrix sense by the diffusion coefficient 
that governs the evolution of the symmetric process and that it depends smoothly 
on the density |12| . 

In contrast with interpretations (a) and (b), the third one is meaningful in di- 
mension d > 2. It has been proved in £Q for asymmetric zero range processes. 
The purpose of this paper is to give a rigorous proof of the third interpretation in 
the case of asymmetric exclusion processes in dimension d > 3. The proof in this 
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context is much more demanding because the process is nongradient. In particular, 
we obtain a non-trivial diffusion coefficient. 

2. Notation and Results 

Fix a finite range probability measure p(-) on 7L d . The exclusion process evolving 
on the discrete torus T d N — {0, . . . ,N— l} d associated to p(-) is the Markov process 
on the state space Xjy — {0, 1} T " whose generator Ln acts on a local function / 

as 

(L N f)(ri) = P(y)v(x){l-v(x + y)}{f(a x - x+ y V )-f( V )} , (2.1) 

x,y£T d N 

where a x,x+y i] is the configuration obtained from r\ by exchanging the occupation 
variables rj(x), rj{x + y): 

{r](z) if z 7^ x, x + y , 

r)(x) if z = x + y , 

r](x + y) if z = x . 

Fix a in (0, 1) and denote by v~ the Bernoulli product measure on Xm with 
density a. Let L* N be the adjoint of Ln in L 2 {v^). This operator is obtained by 
replacing p(y) by p*(y) = p(-y) in (JOJ. 

Denote by T d the d-dimensional torus. Fix a continuous function po : T d [0, 1] 
and denote by v^ r.\ the product measure on {0, 1} T « associated to p n . This is the 

Bernoulli product measure on {0, 1} T ™ with marginals given by 

v po{-M x ) = = Po(x/N) 

for x in T^. 

For N > 1 and a configuration ry, denote by K N (rf) the empirical measure asso- 
ciated to 77. This is the measure on T d obtained by assigning mass N~ d to each 
particle of 77: 

n N (r,) = N- d r)(x)6 x/N , 

x&T d N 

where 5 U stands for the Dirac measure on u. It has been proved in that if 
particles are initially distributed according to v p Q t.\ for some profile po : T d — > [0, 1], 
then tt n (ijtN) converges in probability to p(t, u)du, where p is the entropy solution 
of the Burgers equation 

d t p + q ■ VF{p) = , (2.2) 

where F(a) = a(l — a) and q £ M. d is the mean drift of each particle: q = z p( z )- 
In this article, we investigate the diffusive behavior of the empirical measure ir N , 

that is, its evolution in times of order N 2 . 

As time increases, the solution of Burgers equation (|1.2fl converges to a stationary 

profile which is constant along the drift direction: 

lim p(t,u) = poo(u) = / p Q (u + rq)dr , 
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provided po stands for the initial data. The limit should be understood pointwisely. 
In particular, in a time scale of order TV 2 , the profile of the empirical measure should 
immediately become constant along the drift direction. 
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We shall therefore assume that the initial state is a product measure v , > asso- 
ciated to a profile po constant along the drift direction: 

q-V Po (u) = (2.3) 

for all u in T d . Assume furthermore that the profile is bounded away from and 
1: 

So < po(«) < I -So (2.4) 

for some So > 0. 

Theorem 2.1. Assume that the initial state is distributed according to v p i.\> where 
the profile p satisfies \2.ty . {2.$ . There exists a smooth matrix-valued function 
a(a) = {aij(a), 1 < i,j < d} with the following property. For eacht > 0, tt n (i] t N 2 ) 
converges in probability to p(t, u)du, where the density p is the solution of the par- 
abolic equation 

d tP = Yji,j d *i{ a i,j{p) d VjP) > (n 

p(o,0 = po(0- { ' 

In this theorem, aij(a) = Dij(a) + (l/2)(l — 2a)ai t j, where Dij(a) is the matrix 
given by (|5.9(l and a^j the covariance matrix of the transition probability p(-): 



2/i Vj ■ 

Notice that by the maximum principle, So < p(t, u) < 1 — So for all (t, u). Moreover, 
the solution of the hydrodynamic equation is constant along the drift direction, 

d 

^2li{d Ui p){t,u) = 

because so is the initial data. 

This theorem is an elementary consequence of the following estimate on the 
relative entropy of the state of the process with respect to a local Gibbs state. For 
two measures p,, v on {0, 1} Tjv , denote by H^{p\v) the relative entropy of p, with 
respect to v. 

H N {n\v) = sup | J f dp - log J e f dv} , 

where the supremum is carried over all bounded, continuous functions, which in our 
finite setting coincide with all functions. For t > 0, denote by Sj^ the semigroup 
associated to the Markov process with generator (|2.1|l speeded up by TV 2 . 



Theorem 2.2. Under the assumptions of Theorem \2.1\ on the initial profile po, let 
{pN, N > 1} be a sequence of probability measures on {0, 1} T « whose entropy with 
respect to v^, -, is of order o(N d ): 

HnM^.)) = o(N d ). 

Then, for every t > 0, the relative entropy of the state of the process at time tN 2 
with respect to v^u \ is also of order o(N d ): 

H^SFl^t,.)) = o(N d ), 
provided p(t, u) is the solution of i2.5\) . 
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In view of this result, we can weaken the assumptions of Theorem l2.1l and assume 
ly that the initial state has relative entropy of order o(N d ) with respect to ^(.)- 

3. Relative entropy estimates 



We introduce in this section some auxiliary measures which will play a central 
role in the proof of Theorem 12.21 The statements presented here appeared essen- 
tially in the same form in 4 and We include their proof in sake of completeness. 

Fix a profile po constant along the drift direction and bounded away from and 
1 as in (|2.4() . Denote by p(t, u) the smooth solution of the parabolic equation l|2.5|l . 
Fix < a < 1. For N > 1, denote by f t N the density of fj, N S? with respect to u%. 
An elementary computation shows that is the solution of 

dtf? = N*L* N f t N , 

where L* N is the adjoint of the generator Lm in L 2 (v^). 

Denote by 5" the space of functions f : [0, 1] x {0, l} zd — > R such that 

(1) There exists a finite set A such that for each j3 in [0, 1] the support of f(/3, •) 
is contained in A. 

(2) For each configuration r), f(-, 77) is a smooth function. 

(3) For each density /3, the cylinder functions f(/3, •), fi(/3, ■) have zero mean 
with respect to v$. Here, fi (/?,-) stands for the derivative of f(j3,r)) with 
respect to the first coordinate. 

Let A : R + xT^M be defined by 

, , , , pit. u)(l — a) 

x(t, U ) = tog J \, ,i ■ 

a[l — p{t, u)\ 

A(t, u) is well defined because the solution p(t, u) of the hydrodynamic equation 
(|2.5|l is bounded away from and 1. Denote by ip^ 1 (77) the density of v^, t % with 
respect to v a : 



( V ) = ±- cxp{ J2 Ht,x/N)v(x)} , 



where Z t is a renormalizing constant. 

For functions f; in 3", 1 < i < d, a time t > and integers £ M <C N, define 
the density ij)^{rf) = ip t J'(r)) with respect to the reference measure i/j^ by 

<ffa) = ^T cx p{ E Mt,x/N) V (x)} x 

d 

x exp{- A- 1 ^ Yl duAt,x/N)— £ f,(^ M (x), r x+y r,)] , 
<— 1 x£f d N ye\> 

where Z\ is a renormalizing constant, A a- = {— K, ...,K} d is a cube of length 
2AT + 1 centered at the origin, rj K (x) is the mean density of particles in 1+ A^-: 

and £' = I — A for a finite constant A chosen for the support of UiP^yV) t° be 
contained in Af for all 1 < i < d, \y\ < £ — A. Throughout this article, A stands 



6 



C. LANDIM, M. SUED AND G. VALLE 



for a finite integer related to the support of the transition probability p(-) or to the 
support of some local function. 

In the following, we will need to take M as a function of N and £ as an indepen- 
dent integer which increases to oo after N. In fact we will require M to be such 
that 

Um 1M = , lim — ^— = . (3.2) 

JV->oo N N~>oo M\A M \ y ' 

We present three elementary results which illustrate some properties of the den- 
sity ipt^fiv)- Denote by Sf the smallest integer m with the property that the common 
support of the local functions fi(j3, •), 1 < i < d, < j3 < 1, is contained in A m . 

Lemma 3.1. Assume that Sf < £ < M and that limAr^oo \Am\/N = 0. Fix a 
density f with respect to the reference measure . There exists a finite constant 
C, depending only on f and p(t,u), such that 

,N \ ir__rt\j.N\\ ^ rn\jd-l 



!<(*,.)) " H N {f\^)\ < CN a 



for all N > 1. 



In the statement of this result and frequently in this article, if measures (i, v have 
density /, g with respect to the reference measure , to keep notation simple, we 
denote by H^(f\g) the entropy of f dv^ with respect to g dv^ and by Et[-] the 
expectation with respect to f dv^ . 



Proof. Fix a density /. By the explicit formula for the entropy, the difference 
h nU I - H N (f | is equal to 

ib N r 7 f 



|/log|^ = O(N^) - log | • 

In particular, we just need to show that the second term on the right hand side is 
absolutely bounded by CN^ 1 . By definition of the renormalizing constant Z\, Zt, 
the logarithm is equal to 

d 

logE v M t>) [exp{ - N- 1 ]T d Ul \(t,x/N)(A e U)(v M (^),r xV )}] , (3.3) 

»=1 xGT d N 

where, for a function f in $ and a positive integer £, 



1 



By Jensen inequality, l|3.3(l is bounded below by — CN' 1 ^ 1 . On the other hand, 
since £ < M, Affi(?7 M (0), 77) depends on the configuration 77 only through n(z) for z 
in Am- In particular, since v^ t s is a product measure, by Holder inequality, (|3.3|) 
is bounded above by 



d|A M | f-( ' 



exp 



{ - d\ A M \N- 1 (G> Ui A) (i, z/iV) (A e f< ) (ry M (x) , t xV ) } 



Since by assumption |Am|-ZV -1 vanishes as N j 00, we may expand the exponential 
up to the second order to show that this expression is less than or equal to CN d ^ 1 . 
This concludes the proof of the lemma. □ 
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Taking / = ipK in Lemma 13.11 we obtain a bound on the entropy of r/^Y i with 



respect to 1/ 



AT 
p(«, 



Corollary 3.2. Under the assumptions of Lemma \3.1\ there exists a finite constant 
C, depending only on f and p{t,u), such that 



^«fl^t,.)) ^ CN 



d-l 



for all N>1. 



Corollary 3.3. Fix a smooth function H : T d — > R and a function g in There 
exists a finite constant Co, depending only on f, q, H and p(t, u), such that 



E. 



N- d J2 H{xlN){Aa)(rF{x),T x n) 



E,,i 



N- d H(x/N)(Ai S )( V M (x),T x r))] | < C ^\A M \/N . 



Proof. By the entropy inequality 



E 



N- d H{x/N){A t Q)( V M (x),T x n) 



is less than or equal to 

H „ t j i log E..n 



exp{ 7 iV- 1 Y, H(x/N)(A m )(r, M (x),T x r,)} 



for every 7 > 0. By Corollarv l3.2l the first term is bounded by Oy -1 . On the other 
hand, repeating the argument presented in the proof of Lemma 13. II we show that 
the second term is less than or equal to 



E„ 



N- d Y H(x/N)(A eB )(n M (x),T xV ) 



Cj\Am\ 

_C/,.JV 

TV v Kt.) 



N- d H{x/N)\Ata){jli M {x),T x f,f 



provided that 7|Am|-/V _1 vanishes as N | 00. In this formula, C is a finite constant 
which depends on g and H. In particular, the difference appearing inside the 
absolute value in the statement of the corollary is less than or equal to 



C C 7 |A M | 

7 



N 



Taking 7 = ^N/\Am\, we show that this expression is bounded by C^/\Km\]N. 
Replacing H by —H, we we conclude the proof of the corollary. □ 
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4. Proof of Theorem 12.21 

We prove in this section Theorem 12.21 In view of Lemma 13.11 Theorem 12.21 is a 
consequence of the following result. 

Proposition 4.1. Fix a measure fi N such that Hn(/j n \Vp r.\) — o(N d ). Assume 
that the profile po satisfies 12. 'ip . \2.4P - There exist sequences {fi, n ,« > 1}, 1 < i < 
d, of functions in $ such that 

lim lim sup lim sup N~ d HN (fi N | ip^t ) = 

n-K» tv^oo 

for every t > 0. In this formula, f n = (fi, rl , . . . , fd,n)- 

The proof of Proposition ^. H is divided in several steps. To keep notation simple, 
denote by H' N (t) the relative entropy of fi N with respect to iftK dv^ ■ 

In view of Lemma 13.11 and of Gronwall inequality it is enough to show that for 
every t > 0, 

Hl(t) < o(iV d ,f) + 7- 1 fdsH^S^u"^) (4.1) 

Jo 

for some 7 > 0. Here, o(N d , f) stands for a finite constant such that 
lim limsuplimsupiV~ !i o(A rd ,f„) = 0. 

n^oc N^oc 

The sequence {fi, n , n > 1} is given by Theorem 15. II To keep notation simple, 
we perform all computations for a single function f = (fi, • • • jfd) and then replace 
it by the sequence f n . 

Recall that M depends on N through the relations 1|3.2|) and that I is an integer 
independent of TV which increases to infinity after N. To prove l|4.1J) . we start 
computing the time derivative of the entropy H^ N (t) . On the one hand, a celebrated 
estimate of ^3] gives that 

-Hl(t) < J f» - 9 t log« f )} du£ . (4.2) 



d_ 

It 



On the other hand, a straightforward computation, presented in section |SJ shows 
that the expression inside braces in the previous integral is equal to 

d 

^E E (d Ui X)(t,x/N){r x w; - L* n (AzU)(v M (x),t x v)} (4.3) 

i=l x£T d N 

d 

+ (V2) E E (dl, Uj X)(t,x/N)r x G t Av) 

+ (1/2) E E (^bxMV^fozMTxHijto) 

- ^ (9tA)(t, x/N)r)(x) + J2 (dt\)(t,x/N) v (x)} + o(N d ) . 

x£T d M xGT d M 



ASYMMETRIC EXCLUSION PROCESSES UNDER DIFFUSIVE SCALING 



9 



In this formula, o(N d ) is a term of order N d iM~ 1 < N d , E^n stands for the 
expectation with respect to ip^. dv^ , W* is the current in the «-th direction for the 
adjoint process and Gij(rj), Hij(j]) are local functions given by: 

W* = X>*(y)i/i»7(0)[l-»7(iO], Gijfr) = £>*(y)tfiyi»?(0)[l -»?(»)] > 

y£Z d y<=Z d 



H ij(v) = 2^ p*(y)v(o)[ 1 - v(y)]{yi - Vo,y T n(v M (x) . )} x 



Here and below, is the operator defined by 

(Vx„/)(>7) = f(o*' y v) ~ M 
and, for a local function /i, T/j is the formal sum 

T/> = ^ r x h . 

Since h is a local function, even if the sum of translations is not defined, the gradient 
Wo iV Th makes sense because only a finite number of terms do not vanish. 

We consider separately the sums in 1)4.3(1 . The goal is to replace each one by 
a simpler expression and a remainder denoted by o(N d ). The remainder o(N d ) 
stands for an expression which may depend on time and on the configuration but 
such that 

rt 



lim limsupA^ / dso{N d )jfdv^ 

t^oo N^oo Jo 







for every t > 0. If the remainder vanishes only after taking the limit in f„, we 
denote it by o(f n ,N d ) and we require 



lim lim sup lim sup N 

n^oc at^qo 



-d 



dso(f n ,N d )f?dv» 







for every t > 0. 

We start with the last term of (|4.3() . By Corollary 13.31 we may replace the 
expectation with respect to ip^ dv^ with an expectation with respect to v^ t ^ 

paying a price of order JV°i/|Am |/JV. After this modification, the last line of 14.3fl . 
becomes 

- Y, (d t X)(t,x/N){rj(x)-p(t,x/N)} + o(N d ) . 

x£T d N 

Since d t \ is a smooth function, we may further replace r](x) by r) {x) paying a price 
absolutely bounded by C£ 2 N d ~ 2 for some finite constant C. 

To estimate the order N d+1 term of (|4.3|l . we first take advantage of the assump- 
tion that the solution pit, u) is constant along the drift direction. 

By paying a price of order 0(t 2 N d ^ 1 ), we may replace the current W* by an 
average \Ag | _1 J2 y eA e T v^i ■ Here again one should keep in mind, that the average 
is in fact carried over a cube of length slightly smaller than 2£ + 1 to ensure that 
all local functions T y W* have support contained in Ag. 



10 



C. LANDIM, M. SUED AND G. VALLE 



Recall that q — (gi, • • • , <7d) denotes the drift of particles. The average of the 
current W* can be written as 

-^Ti E T ' W i = ?*{/(0)-2r/(0)/(0) + r/(0) 2 } + E <(^(0),r,»?) , 

where g* = —qt and 

= - E fo(°) ~ a l t 7 ?^) - «] - « E P*(f)f< iv(y) - ■ 

y£Z d 



The first term of the current gives no contribution since for any function J, 

d 

^E E (duMt,x/N)qiJ(/(0W(x)) = 

*=1 xeT d N 

because Ei<i<d9i( 9 «*^)(*)") = {p( f > «)[! ~ Ekk^ ftftirft*. «) van- 

ishes for all (t,u). The first term of l|4.3|) becomes therefore 

d 

^E E iduMt,x/N)r x {(A e w*)(v\0),v) - ^(Aif^iOU)} . 

To ensure that the function which appears in Aiw* has mean zero with respect to 
the all canonical measures on the cube Ai, we further replace Agw* by A eW* , where 

{A\w*){a,n) = {A e w*)(a, v ) + Qi ^ ~ 7 ■ 

This replacement is pemitted because qid Ui p = 0. 

Following the nongradient method, we add and subtract Xa<j<d Dij(if N (0)) 
[if N { e j) — ^^(O)]- Since the diffusion coefficient is smooth, this expression is equal 
to ^KjKdidijiv 61 * 1 ( e j))~ <i,i07 EJV (0))} + ( £ Af) _2 > where dij stands for the integral 
of Dij. In particular, after a summation by parts, the first line of (|4.3|) . may be 
rewritten as 

d 

N J2 E (duMt,x/N)T x Vf N ' M < e ( V ) (4.4) 

d 

where 

d 

(A° lW *)( v e (0), v ) + J2 D ^(v K (0))[v K (e 3 )-V K (0)} L* N (AtU)(ri M (0),v)- 
3=1 

It is not difficult to see that there exists a finite constant C(a) such that 
H (fi N I v^) < C(a)N d for every probability measure fi N on {0, 1} T «. In par- 
ticular, by the usual two blocks estimate, since dij is Lipschitz continuous, for 
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every T > 0, 

lim lim sup lim sup 

f dt [ ^(d V )f t N (v)N- d \ d % ,M N {x))~d t ^{x))\ = 0. 

We may therefore replace in the second line of (|4.4|1 the average of particles over a 
small macroscopic cube by the average over a large microscopic cube, i.e., replace 
rf N (x) by T] £ (x). 

On the other hand, the usual nongradicnt techniques, based on integration by 
parts formula, allows the replacement in 1)4. 4[) of Dij(r] eN (0))[rj £N (ej) — rj eN (0)] by 
Aj(^(0))[^ (e 3 ) - rf'(p)]. Here I' = I - 1 for the previous function to depend 
only on the sites in An. To keep notation simple, we will denote this expression by 
Dij(v e (0))b/( e j) — ^/(O)]. We refer to Chap. 7 of |H| for a proof of this replacement. 

In subsection !6 . 21 we prove that we may replace L* N (Aifi)(r] M (0), rj) by L* A (Atfi) 
(»/ (0), 77). Here L\ stands for the restriction of the generator L* N to the cube A^. 
This means that we suppress all jumps from Ag to A£ and all jumps from A^ to 
Ai. In particular, this generator leaves rj e (Q) invariant and it is acting in fact only 
on the second coordinate. This replacement is one of the main technical point of 
the article. It is here that the special form of ipK plays an important role, that we 
need the spatial averages and the particular size of M and £ presented in (|3.2(l . 

Up to this point, we transformed the first line of l|4.3|l in 

d 

N^2^(d u Mt^/N)r x Vf( V ) (4.5) 

d 

+ E E (^axwwjOA*)) + o(N d ) , 

where 

d 

Vf( V ) = (^<)(^(0),rj) + E^(^(°))[ ? Ae,)-^(0)] - L^f 4 )(?/(0), 77) . 

i=i 

By the nongradient method, the first line can be shown to be of order o(f, N d ). 
Details are given in subsection 16.41 

It remains to consider the second and third line of l|4.3l) . By the one block 
estimate the second line of (I4.3H is equal to 

d 

(1/2) £ £ {d^. U] X)(t, x/N) Oi^T x F(rj l (0)) + o(N d ) , 

i,j=l x£T d N 

where tr^j is the symmetric matrix defined just after (|2.5|) and F(a) — a(l — a). 
For 1 < i, j < d, let 



Jij(fi) = 2/3(1 - 0){Dv(P) ~ faij} 



(4.6) 
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We prove in subsection 16.31 that the third line of 1(4.3(1 is equal to 

d 

(V2) E E {d^\){t,x/N)(!^\)(t,x/N) ai jFtf{x)) 

d 

+ (V2) E E (^bxMfaiV&xMJutfi*)) + odn,N d ) . 

In conclusion, we proved that (|4.3|l is equal to 

2 d 

E E E G™.(t,^/iV)ff™.(^(x)) (4.7) 

m =l i,j=l xgi^ 

- E {Bt\){t,x/N){rf{x)- p{t,x/N)} + o(f n ,N d ) . 

where 

Gjj(t,u) = (aS liUj A)(t,tt) , fl^CS) = d^OS) + (l/2)a iJ F(/3) , 

Gf^t.u) = (9„ 4 A) (t, u) (d Uj A) (t, u) , Hf d {p) = (l/2){J h3 (p) + a^F(f3)} . 
An integration by parts shows that 

2d 

E E / G™ (t, u)H™ (p(t, u)) = 0. 

In particular, in formula 1)4.7(1 . we may replace the terms ff™ (r/ f (x)) by iJ™(?/(:r))— 
H™(p(t,x/N)) paying a price of order o(N d ). A further elementary computation 
gives that 

2 <Z 

E E G&(*,tt)TO)'(p(*,«)) = ($A)(t,«) 

m— 1 — 1 

for every i and u, where (H^)' stands for the derivative of ffjj}. Therefore, (|4.7() 
becomes 

2 d 

E E E G^(t,x/N)B^( v £ (x),p(t,x/N)) + o(f n ,N d ) , 

m=l i,j=l igTj, 

where 

5™ (a, 6) = - flJ3(6) - )'(b)[a-b] . 

At this point we may repeat the standard arguments of the relative entropy method 
do conclude. We refer to Chap. 6 of for details. 

5. HlLBERT SPACE OF VARIANCES 

We prove in this section the existence of functions f i , . . . fd in J which approxi- 
mate the current in the Hilbert space of variances. We rely on recent results based 
on general duality presented in ^D], |12|. 

For < a < 1, denote by Q a the space of cylinder functions g such that E Va [g] = 
3 a E Va [g] =0: 



g(a) = E Va [g] = and g'(a) = -j^E Vf3 [g\ 
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For each function g in Q a we define by 

Nil* = \g\l + \\g\\\ a , (5-1) 



where 



g\l = sup {2 YV V x, < g: r](x) > a -^^a ■ era) 
„ = sup {2 < a, h > Q - < h, {-L s )h > Q } . 



\g\\\ 

In this formula, x( a ) = a (l — a )i a ' b stands for the inner product in M. d and 
<C •, • 3> Q for the inner product in Q a given by 

< g,h > Q = ^ < g; r x h> a , 

where < f\ ; fa > a denotes the covariance of /1, f% with respect to v a . Notice that 
in the sums which appear in the formulas above, all but a finite number of terms 
vanish because u a is a product measure. Theorem 15.11 is the main result of this 
section. 

Theorem 5.1. There exist a smooth matrix-valued function D(a) = {Dij(a), 1 < 
h j < d} and a sequence of functions {fi, n , n > 1} in 3, 1 < i < d, such that 

d 

lira sup flu;? (a, rj) + } Dij(a)[r)(ej) - rj(0)] - i*f i n (a,7j)|| Q = 



n ^°°ae[o,i] 



/or 1 < i < rf. Moreover, for any vector v in M. d , 

d d 

lim Y] < y2v j f j , n (a,r)),(-L s )T x y2v j f j , n (a,ri) > a (5.2) 

n — >oo L — * L — * L — * 

xei d 3=1 3=1 
= xi a ) v ' {D{a) — aa}v 



uniformly in a. 

This result is a slight generalization of Corollary 10.1 and Lemma 10.4 in [5], 
proved in |5J using results presented in We have the advantage here to obtain 
uniformity up to the boundary. In sake of completeness, we present a simpler proof 
relying on the generalized duality developed in ^U], |12) . 

To keep notatiom simple, we prove Theorem 15 . II for the current Wi obtained from 
w* by replacing p*(-) by p(-) and for the generator L in place of L* . 
Duality. For each n > 0, denote by £ n the subsets of Z d with n points and let 
£ = U n >o£ n be the class of finite subsets of I/ 1 . For each A in £, let ^a be the 
local function 

By convention, — 1. It is easy to check that {"P a, A S £} is an orthonormal 
basis of L 2 (v a ). For each n > 0, denote by V n the subspace of L 2 {v a ) generated by 
{$>A, A £ £ n }, so that L 2 {u a ) = (B n >o'Dn- Functions in V n are said to have degree 
11. 
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Consider a local function /. Since {^a ■ A € £} is a basis of L 2 (v a ), we may 
write 

/=EE k«»^)*a. 

n>0 AGfr, 

Note that the coefficients f (a, A) depend not only on / but also on the density a. 
Since / is a local function, f : £ — > M is a function of finite support. 

Fix a local function / and denote by f(a, A) its Fourier coefficients. / has zero 
mean with respect to v a if and only if f(a, cf>) — 0. It belongs to Q a if and only if 
f (a, <fi) = and the degree one part is such that 

f(a,to) = 0. 

In this case, we may rewrite the degree one piece as 

V^^ftMz})^)-^)]. 

In conclusion, all functions / in Q a may be written as 

^E^wiiiW-ii + E E K«.^)*A. 

zGZ d n>2Ae£„ 

For n > 0, denote by 7r„ the projection on T> n so that / = X) n >i 7r «/ f° r / m 
Q a . In the formula above, the first term corresponds to nif, the piece of / which 
has degree one, and the second term corresponds to (7 — 7Ti)/, the piece of degree 
greater or equal to 2. 

It is clear that a local function of type h — r x h belongs to the kernel of the 
inner product <C •, • ^> a defined above. This is the case of rj(z) — 77(0) so that 
|| /|| -i, <* = ||(7 — 7Ti)/||-i,a- In contrast, any function h of degree greater or equal 
to 2 is such that 

y~]xj<h; T](x) > a = 
for all i so that \h\ a — 0. Therefore, |/| Q = |7ri/| Q and 

for every local function / in Q a . 

The generator on the Fourier coefficient. Let £* be the class of all finite 
subsets of = Z d \{0} and let £„^ n be the class of all subsets of with n points. 
For a local function / in Q a , define If : [0, 1] x £» — > R by 

(%f)(a,A) = Ef(",[AU{0}]+z), 

zez d 

where f(a, B) stands for the Fourier coefficients of /. In this context, a function 
f(a, rf) belongs to Q a if and only if f(a, <f) = (1f)(a, 4>) = 0. It has been in proved 
in |12| that for every zero- mean local functions /, g 

« f,g »a = < (If), (?9) > - E — T7 E WX«^) (^9)(a,A) . (5.3) 

* — ' n + 1 z — ' 

n>0 AGS,, n 

For functions in Q a , this sum starts from 1 because (%f)(a,<p) — (%g)(a, cf>) = 0. 

Observe that not every function f : [0, 1] x £« — > ffi is the image by 1 of some 
local function / since 

(T/)(a,A) = (Tf)(a,S z A) (5.4) 
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for all z in A. Here, S Z A is the set defined by 

_ / A-z iizgA, 
2 ~ 1 {A-z)o,-. HzeA. 

Let f* : [0, 1] x ^ 1 be a finitely supported function satisfying l|5.4|l . Define 
f : [0, 1] x £ -> R by 

f(a , B) _ { V{ o» .^o, m 

An elementary computations shows that 1f(a,r)) = f*, if f(a,rj) is the local func- 
tion whose Fourier coefficients are f(a, A). Notice that f(a,rj) belongs to Q a if 
f») = 0. 

For any local function /, %(Lf) = £ a T/, provided 



£ Q = £ s + (1 - 2a)£ d + + £-} 

and, for A e £*, t> : £* — > M a finitely supported function, 

(£.0)(B) = (1/2) °(V -9(B)] + J2 s (yMSyB)-9(B)], 

x,y&i y&B 

(£ d 9)(A) = a(y -x){9(A x , y )- 9(A)) + J2<y)MS y A)- 9(A)}, 

x£A,y£A y&A 
x,y^0 y^O 

(£+9)(A) = 2 a(y-x)9(A\{y}) 

x£A,y£A 

+ 2 ]T a(x){9(A\{x}) - 9(S x [A\{x}})} , 

x£A 

(£_t))(A) - 2 £ a(|/-x)t)(4uM). 
In this formula, A^y is the set defined by 



(A\{x})U{y} ifxeA,y?A, 
A x , y = <( (A\{y}) U{x} if y e A, x $ A, 
^4 otherwise ; . 

Hilbert spaces. For two local functions /, g, let 

< /, 9 >a,i = < /, (-i s )9 >a 

and let -ffi(a) be the Hilbert space generated by local functions / and the inner 
product <C •, • S>q.i- Denote by <C •, ■ 3>i the scalar product on defined by 

<f,0>i=y"^-r $2 f(a,A)(-£ s0 )(a,A) 
n>0 Ae£,, n 

and by S)\ the Hilbert space generated by the finite supported functions endowed 
with the previous scalar product. From the previous definitions, for every local 
function /, g, 
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To introduce the dual Hilbert spaces of Hi, Sji, for a local function /, consider 
the semi-norm || • ||_i given by 

ll/ll-i,a = supf2< /,g> Q -<g,g^>i, a } , 

where the supremum is carried over all local functions g. Denote by the Hilbert 
space generated by the local functions and the semi-norm || • \\-\. In the same way, 
for a finitely supported function f : £ * — ► R, let 

Hfllij = sup{2 <f,fl> - < 0,g >i } , 

where the supremum is carried over all finitely supported functions g : £ * — > K and 
< ■, • > is the inner product on L 2 (£*) defined in 1|5.3|) . Denote by S)-\ the Hilbert 
space induced by the finitely supported functions f : £ * — ► K and the semi-norm 
|| • || _i. By the identities for the L 2 and the Hi norms, we obtain that 

\\f(<X,V)\\-l,a = II(£/)(<V)H 2 -1 (5.6) 

The currents. Recall the definition of the current Wi(a, rf) given in sectional Wi 
is obtained from w% by replacing p*(-) by p(-) and can be expressed as 

w t = -a(l - a) P(y) Vi *o,y ~ " X! p ^ ^ ^{y) ~ v(fy} ■ 

Denote the first piece, which has degree 2, by a(l — a)w®. On the other hand, since 
v( e k) — ^(0) = r)(ek + x) — i](x) for the norm | • | Q for any x, the piece which has 
degree one is equal to aX) yeZ<i J2i<]<dP(v) yiVji^j) ~ v( )} so tnat 

d 

Wi = a(l - a)w° - tt^gy [y(ej) - t?(0)] . 
i=i 

Let fdi = An straightforward computation gives that 

tDi(a,{z}) = -2z, ; a(z) 

for z ^ and tt>i(a,A) = otherwise. Notice that XVi does not depend on a. We 
have now all elements to prove the main result of this section. 
Proof of Theorem 15.11 Fix 1 < i < d. By Theorem 4.1 in 12 , ro; belongs to 
f)_i because Wj(a, (p) = and we are in d > 3. 

It has been proved in Lemma 4.3 of ^21 that for each A > there exists a solution 
fi.\(a,A) of the resolvent equation 

Afj,A — £afi,\ = tt>i 

satisfying (|5.4|l and such that fi,\(a, 4>) = 0. 

By Theorem 4.4 in \T2\ , for any k > 1, there exists a finite constant CV- indepen- 
dent of a and A such that 

A^(l + n) fe < 7r n fi,A,7rnfi,A >„ + ^(l + n) fe < 7r„fi,A, (-£ a )7r n fi,A >„ < C fc 

n>0 n>0 

. (57) 

for every A > and a in [0,1]. In this formula, n n stands for the projection on 
£*,n- (ftnf) A) — f(a, A)1{A € £*,n}i an d < •, • > n for the inner product in 
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with respect to the counting measure: 

< f,fl >n = K a ,A) Q (a,A) . 

The estimate is uniform in a because the current fo^ does not depend on a. 

By section 6 of f° r each z in Z d , fi,\(-, {z}) is a smooth function in [0, 1] 
and there exists a subsequence Xk | such that U,Xk( a i{ 2 }) converges uniformly, 
as well as its derivatives, to some smooth function fi(a, {z}). 

By the proof of Lemma 2.8 of taking a further subsequence, we may assume 
that — (£ a fi,\ k )(a, •) converges weakly to Wi in Sj-i for a countable dense subset of 
densities a in [0, 1]. By Lemma below. — (£ a U,x k )( a : converges weakly to to.; 
in Sj-i for all a in [0, 1]. 

Our goal is to replace the sequence fi \ h by a sequence t)i^ n of finite supported 
functions with all the above properties of fi t \ k and for which — (£ Q fji jTl )(a, •) con- 
verges strongly to Wi in uniformly in a. 

For each a fixed, we may take convex combinations of the functions fa \ h to 
obtain a new sequence Qi^ n such that — £ a 8i,n converges strongly to Wi in 
Lemma 15 . 21 below shows that the procedure can be made uniform in a. Indeed, fix 
e > and a finite set {aj, 1 < j < m} in [0, 1]. The standard procedure to derive a 
strong converging sequence from a weak, bounded converging sequence shows that 
there exist M > 1 and a probability (0%, . . . , 9m) in {1> • • • , M}, such that 

max ||£ a ,0(0,-, •) + w t \\-! < e, 

l<j'<m 

where 

M 

fl(«J.-) = 53^fi,A,(Q!j,-) ■ 
1=1 

Notice that we are taking the same convex combination for all densities aj . If m is 
equal to <5 _1 , given by Lemma f5 . 21 below . and aj — jS, by Lemma f5. 21 

sup ||£ a fl(a, •) + K>t II -1 < 2 £ , 

a£[0,l] 

where g(a, •) is obtained from fi t \(a, ■) through the same convex combination. We 
have thus constructed a convex combination which guarantees the strong conver- 
gence in for all values of a. That is, there exists a sequence Qi >n { a i') such 
that 

• For each n > 1, and each z in Z d , 0i, n ('i{ z }) is a smooth function of a 
which converges uniformly, as well as all its derivatives, to some smooth 
function fi(a, {z}). 

• Each Qi :1l satisfies (15.411 and Qi yn (a, cf>) — because the functions fi t \ h have 
this property. 

• The sequence converges uniformly to Wi in 55 _i: 

lim sup \\£ a Bi,n(a, ■) + tt)i||-i = 0. 

It remains to replace the functions by finite supported functions. Fix two 
integer m, I and let E}j lTl (a!,A) = Qi, n (ot, A)1{\A\ < m,A < A(}. The integers m, 
f , which depend on n and increase to infinity with n, will be chosen later. Here, 
A < if there exists z in A such that 5* 2 A C A(. In this way, f)^ satisfies (15.411 . 



18 



C. LANDIM, M. SUED AND G. VALLE 



The sequence t)i >n just defined has the first two properties of the sequence Qi, n 
enumerated above because m and £ increase to infinity as n f oo. To prove the 
third one, recall from the computations performed after (4.12) in |§] that 

||£afli,n(a, •) - £ a t)i,n(®, II —X 
m+1 

< C ^2(l + k) \\n k Q hn (a, •) - 7T k hi,n( a > Ollo.fc 
fc=i 



^(1 + A;) 2 < 7r fe n i: „(a, •), (£ s 7Tfc0i,n)(a, •) > fc 



for some finite constant Cq independent of a. Here, || • || o,fe stands for the norm 
associated to the scalar product <•,•>& defined above. By (|5.7|l . the second term 
on right hand side can be made uniformly small in a by choosing m large enough 
because each function is obtained as convex combinations of the solution of 
the resolvent equation. For a fixed finite set ot\, . . . ,a r , we may turn the first term 
as small as one wishes for {a^ 1 < i < r} by taking I large enough. By Lemma 
15.21 below, we may turn the estimate uniform in a because the functions Qi, n are 
convex combinations of the solution of the resolvent equation. 

For each fixed n, the functions f)i.„ (a, •) has a uniform support. Since \)i in satisfies 
(|5.4(l and f)i ira (a, <j>) = 0, the local functions fi, n {&,v) obtained from f)i.„ through 
(|5.5(l are in 

We claim that the sequence — x(&)fi t n( a j v) nas the properties required in the 
statement of the theorem. In view of the decomposition of the current Wi, by Q5.ip . 

d 

Wi(a,i)) +^2,D itj {a)[r){ej) - r?(0)] + x(a)£/ i , 7l (a, r?) 
j=i 

II 2 
= x(") w i ) + x(a)(/-7ri) J L/ ij „(a,r7) 

I — l.a 



(5.8) 



+ |5Z{A,j(Q!)-Q!o-i,j}[»7(ei)-»7(0)]+x(a)7rii/i 1 n(a,»?) ■ 

I — a 

Since functions of degree 1 are in the kernel of the scalar product <C ■ , ■ ^$> a , we 
may replace (I — 7Ti)£/i,n by Lfi^ n on the first term on the right hand side. On the 
other hand, by definition of 1, by identity (|5.t)|) and since Tu>° = K>i, the first term 
on the right hand side of l|5.8|) is equal to 

x{a) 2 \\tt>i + £ Q t)i,n(a, ■)!!_! • 
This expression vanishes, as n | oo, uniformly in a, by construction of the sequence 

On the other hand, an elementary computation, presented just after (5.4) in |12j . 
shows that 

■KiLfi, n (a,ri) = ^2 a(z)t)i yn (a, z)[r](z) - i](0)} . 

z&L d 

Since rj(z) — ^(O) = Xa<j<d z j[v( e j) f° r the norm | • | a , the second expression 

on the right hand side of l|5.8|l is equal to 

d 2 

V{Aj(a) - aa hj + hij, n (a)}[T}(ej) - 77(0)] , 

I * — ' a 
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where 

K,j,n(oi) = x(a) Y a(z) Zjt)i. n (a,{z}) . 

By construction, f)i n (a, {z}) converges to fi(a,{z}) uniformly, as n | oo. In par- 
ticular, if we define Di^(a) as 

Aj(a) = aa itj - x(a) Y a(z) zj f 4 (a, {z}) , (5.9) 

it not difficult to show from the variational formula for the norm | • | Q that the 
second term on the right hand side of (|5.8|l also vanishes as n | oo, uniformly in a. 
This proves the first statement of the theorem. 

Notice that Dij(-) inherits the smoothness of fi(-, {z}). 

It remains to check identity l|5.2[) . By definition of the scalar product <C ■, • ^> a 
and by identity (|5.3I) . for any vector « in R d , the left hand side of i|5.2|) with the 
sequence — x(&)fi,n(&, rj) in place of fi <n (a, rj), is equal to 

d d 
3 = 1 3=1 

Since — £ a l)j t n converges to K>j in uniformly in a, and since f)_y )n is (n, a)- 

uniformly bounded in f)±, the limit of the previous expression is equal to the limit 
of 

d d 

X(a) 2 Y VjV k < t)^„{a,-),n> k > = ~x{a) 2 Y v i Vk X! ^jA a A z }) z ka{z) . 
j,k=l j,k=l zezi 

The last identity follows from the explicit formula for rDfc. Notice that a factor 1/2 
appeared because in the definition of the inner product < •, • >, there is the term 
(n + 1) _1 . By construction, (jj jn (a, {z}) converges uniformly in a to fj(a, {z}). In 
particular, the previous sum converges uniformly to 

d 

-X(a) 2 Y v i Vk Y fj( a A z }) z ka(z) 

j,k=l zEli 

By definition (|5.9|l of the diffusion coefficient Dij(a), the previous expression is 
equal to 

x(a) v ■ {D(a) — aa}v . 
This concludes the proof of the theorem. □ 

We conclude this section with a technical lemma needed above. 
Lemma 5.2. For each e > and k > 1, there exists S > such that 

sup \\£ a fi t x{oi,-) -£/3f;,A (/?,-) II -1 < £ 
\a-/3\<6 

sup A^(l + n) fe Y, {UA^A)-\. hX (f3,A)} 2 < e 

\<*-P\< d „>i Aef», n 

for all < A < 1 . 
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The proof of this lemma is implicit in the proof of the regularity of f\(-,A) 
presented in Theorem 5.1 and Lemma 5.2 of |l()j . We just need to write the equation 
for fi t \(a, •) — fi.\{P, •) as a resolvent equation with a right hand side which is a 
function in Sj-i times 0(a — (3). Details are left to the reader. 

6. Technical bounds 

We present in this section some technical lemmas and some computations omit- 
ted in sectional 

6.1. Computation of N 2 L* N ipj* f /ipj* f . Since L* N is the generator of the exclusion 
process associated to the transition probability p*(y) — p{—y), 

. £ v(*)H-v(* + y)V(y){ " ,» -'}■ 

For each fixed bond (x,y), ipK (<j x,x+v ri)/ip^ (77) is an expression of order iV -1 
because fj(-, 77) is a smooth function for each fixed configuration 77. We may therefore 
expand the exponential up to the second order. The order one term is exactly 
N 2 L* N logoff and is responsible for the first two lines of l|4.3|) plus a remainder of 
order N^ 1 . The second order term is equal to 

(1/2) Y, v{x)[l-v{x + y)]p*{y){N{\{t,x + y/N)-\{%x/N)} 



E E (^A)(M/A0V,,, + ,(A^(^),^)} 



Since I + Sf + A < M, the gradient V x ^ x + y acts either on the first coordinate or on 
the second but never on both. fi(-,rj) being a smooth function, the contribution of 
the gradient V ' x ,x+y applied on the first coordinate is at most of order M~ d . Since 
there are 0(M d ~ 1 ) boundary sites z for which V XtX + y rj M (z) does not vanish, the 
total contribution of the gradient V x , x + y acting on the first coordinate of Aefi is of 
order M -1 . 

We consider now the set of sites z for which the gradient V XiX + y acts on the 
second coordinate of Afc. In this case, z should be at a distance smaller than £ + A 
from x and we may replace (d Ui X)(t, z/N) by (d Ui X)(t, x/N) paying a price of order 
( d+1 N~ 1 . At this point, for a fixed i, after a change of variables z' — z — x,we may 
rewrite the sum appearing inside braces in the previous formula as 

(d Ui \)(t,x/N)T x V ,y ]T r^- UV M (Z\T Z+WV ). 

Since the summation over z takes place on Af +J 4, we may replace r\ M (z) by ?7 M (0) 
paying a price of order £/M. In this case the previous sum becomes 

{d Ui X)(t,x/N)T x W Q ,y E U(v M (Q),t z v) = (d Ui \)(t,x/N)T x Vo,yT U (r, M (o),-) 

because the contribution of each fixed w is the same after replacing r\ M (z) by ?7 M (0). 

To obtain the third line of l|4..'-{(l and the correct order of the remainder, it remains 
to expand N{X(t, x + y/N) — X(t 7 x/N)} and to develop the square. 
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6.2. Replacement of L*(A £ fi)(f7 M (0), 77) by L\ ( (Atji){rf{Q), 77). Observe initially 
that the generator acts either on the first coordinate or on the second but never 
on both because we assumed that sj + I < M. Hence, we have to show that the 
action of the generator on the first coordinate is negligible. This is the content of 
the next result. 

Lemma 6.1. Fix a function f in J, a smooth function G : R+ x T d — > K and 
assume that M satisfies conditions \3.ty) . For every T > 0, 

lim limsup 

' dt f N 1 -* £ G(t,z/N)T z (L*-L* M )(Aef)( V M (0),r))f t N dp»\ = 0. 

Notice that in L A (Aij)(i] M (0), rf), the generator is acting only on the second 
coordinate because £ < M. 

Proof. Let fi(a, rj) = (d a f)(a, rf). Since f(a, •) is a smooth function, the contribution 
of (L* - L* Ae )(A e U)(r] M (0),r!) is equal to 

N i-d M -d J- G(t,z/N)r z v{x)[^-il{x+yW{y){Azh){ri M {Q),ri) + ON(l) 

z£T d N xGA m 

(6.1) 

plus a similar term with a negative sign and x + y in Am, x not in Am- Here 
the remainder ojv(1) is of order N/M d+1 . From this point, the proof is divided in 
several steps. 

Step 1. The first one consists in translating the local functions 77(3:) [1 — rj{x + y)], 
which lies at the boundary of Am, by few steps in order to have their support 
contained in Am- For this purpose, it is enough to show that for every fixed y, 

N i-d M -d G(t,z/N)T, £ r x W(Aih)(V M (fl),v) (6-2) 

x+y£A M 

is negligible if W — h — r ei h for some local function h. Here and below, a function 
HN,e{t,r/) is said to be negligible if 

l-T 



lim limsup 



dt / H NA {t,rf) jf dv* 







for all T > 0. Since there exists a finite constant Co such that Hm{ij. n \ vf£) < 
CoN d for all measure fi N , by the entropy inequality, Feynman-Kac formula and 
the variational formula for the largest eigenvalue of a symmetric operator, to prove 
that a function is negligible, it is enough to show that 

lim limsup / disupj / H{t,rj) f du* - sN 2 - d D N {f)\ < (6.3) 

for every e > 0. Here, the supremum is carried over all densities / and Dpf(f) is 
the Dirichlet form given by Djv(/) =< — Ljvv 7 /, v 7 / >, where < •, • > stands for 
the inner product in L 2 {v^). 

Since the local function W has mean zero with respect to all canonical invariant 
states, W — L A w for some finite set A and some local function w, where L S A stands 
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for the symmetric part of the generator L restricted to the set A. In particular, we 
need only to show that 

N i-d M -d £ G (t,z/N)r z ]T T x CV b w)(Ath)(v M (P),V) 

x+y£A M 

is negligible for a fixed bond b = (pi, 62) and a fixed local function w. Fix < t < T, 
a density / with respect to and consider the linear term in variational formula 
(IQl : 

N i-d M -d J- G(t,z/N) £ J r x {V b w){A^){ n M {Qlri)r. z fdv^ , 

x+y^Au 

where we performed a change of variables £ = t z t]. Since r K Vb ~ ^b+xT x , per- 
forming a change of variables £ = a b+x rj, we may rewrite the previous expression 
as 

N i-d M -d J- G (t, z/N) I T x w(Mi)(v M (0),v)V b+x T- z fd^ 

zeTj x<EA M 
x+y£A M 

plus a term of order NM^ d ^ 1 . This term appears when taking the difference 
V h+a; (^fi)(77 M (0), 77) which is absolutely bounded by CNT d . 

Rewrite the difference a — b — T^ z f(a b rj) — T- Z f(r]) as (^/a — ^/b)(^/a+ y/b) and 
apply the elementary inequality 2ab < -fa 2 + j~ 1 b 2 , which holds for every 7 > to 
estimate the previous expression by Ce~ x M~ 2 + eN 2 ~ d D]y(f). This proves that 
(|6.2(l is negligible, concluding the first step. 

Step 2. Once that all functions have been translated to have its support contained 
in Am, we take advantage of the fact that each function which appears in (|6.1|) 
at one side of the boundary, appears also at the other side with reversed sign. In 
particular, adding the intermediary terms to complete a telescopic sum, after (|6.2|l . 
(|6.1|l can be rewritten as 

m 

i=i zgTj xeA M -A 

for a family of local functions hj — gj — r ei gj for some 1 < i < d. Here m is a 
finite integer which depends on p(-) only. In particular, the local functions hj have 
mean zero with respect to all canonical invariant measures. Here again, A is taken 
large enough for the support of each local function r x hj to be contained in Am- 
We claim that such a term is negligible. 

Since all local functions h which have mean zero with respect to all canonical 
invariant measures can be expressed as L s A h Q for some finite set A and some local 
function /i , fix a bond b, a local function h a and consider the linear term in l|tj.3[l : 

N i-d M -d ^2 G(t,z/N) Y J MVbh )(A e H)(r, M (Q), V )T_ z fd^ . 

z£T% x£A M _ A 
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Since T x ^7b — ^7b+xT x , a change of variables £ = a b+x r], similar to the one performed 
in the first part of the proof, permits to write the previous expression as 
N i-d 



yfil ]T G(M/A0 E J T x h (A e h)(v M (0),« b+x v)V b+x T^ z fd^(6A) 

zeT d „ xeA M -A 

E G(t,z/N) E f^h V b+x (A£h)(v M (0),v)r- z fd^. 



M d 

We claim that both terms can be estimated by eN 2 ~ d Diy(f) and an expression 
which vanishes as N f oo and then I \ oo. Notice that in the second term, the 
gradient Vf,+ X is acting only on the second coordinate. 

Consider the first line of i|6.4[l . Repeating the arguments presented at the end of 
the first step, we may bound this integral by the sum of eN 2 ~ d D^(f) and 

WfisY, E /(^W(^(0),^) a {r-,/W+r_./(^)}^ 

zGTt x£A M -A 



for some finite constant C. Notice that we got an extra factor TV -1 in this passage 
and that we included G and /io in the constant. We perform a change of variables 
£ = (j b+x rj and denote by / the average of the translations of /: / = N~ d J2zeT d Tz -f 
to rewrite the previous sum as 

Ce-'J ^f^ M (0)^) 2 /(r?)<^ + 0(£ d M- d ). 

Here we took advantage of the fact that (A £ fi)(r; M (0), o b+x rj) = (A £ fi)(f7 M (0), 17) 
unless x belongs to Ag. Since f i (- , 77) is a smooth function, uniformly in r\, the 
integral in the previous expression is less than or equal to 

Ce- 1 J {AfaWiUUfmdvZ + Ce- 1 J {^(0) - ^(0)} 2 /(„) di£ . 

The usual proof of the two blocks estimate permits to show that the second integral 
can be estimated by eN 2 - d DM{f) and an expression which vanishes as TV "f 00 and 
then £ | 00. We leave the details to the reader. In contrast, the usual proof of the 
one block estimate permits to show that the limit, as N f 00, of the first integral 
minus eN 2 - d D N (f) is bounded by 

Ce _1 sup / { 1 E T vh(K/\Ae\,v)} d/j, Ae ,K ■ 

In this formula, fJ-Ai.K stands for the canonical measure on concentrated on 
configurations with K particles and the supremum is carried over all integers < 
K < \At\. Divide the average in Ai in two averages and recall from Lemma A. 7 in 
[S] that the Radon-Nikodym derivative d/J,A e ,K/di' K e /\j ii 1 is bounded, uniformly in 

K, provided v^ 1 stands for the grand canonical measure on Ai with density (3. The 
previous expression is thus less than or equal to 

Ce- 1 sup f E T v h(J3,r,)Y <k# . 

0<P<lJ ^iHlly^ j 

In this formula, A^i stands for one half of the cube Ag. Since fi(a, •) is local 
function, with uniform support and which has mean zero with respect to u£j , the 
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previous expression is of order l~ d because i/£ is a product measure. This conclude 
the estimation of the first term in H6.4(l . 

We turn now to the second term of i|tj.4fl . Notice that the gradient Y7 'b+ x {Ag\i) 
(r] M (0), rj) vanishes if x does not belong to Ag + A- In particular, 

T K / lo V 6+;c (A,f 1 )(? 7 M (0),r;) = r 3 ,hoV i+ai (A4h)(Ti M (0) ) T i ) (6.5) 

a:eA M -A x£A e+A 

is bounded by a constant which does not depend on N. On the other hand, for 
every < K < \Am\, repeating the computation presented in the second paragraph 
of the second step, from the end to the beginning, we obtain that 



J2 J ^feoVt+xtAifiK^tO),!))^^ 

J (r x W b h a )(A i h)(v M (0),v)dfiA AuK 



x£Ae +A 



x£A e+A ' 

Summing over all bonds b, we recover L s h,Q = h = g — T ei g, for some local function 
g and some 1 < i < d. The previous expression is thus equal to 



]T J ' (T x g)(A £ U)( v M (0), V )dfi AMtK 

f(r x g)(A e h)(ri M (QU)d f XA M ,K 



where d~Ai + A stands for the lower boundary in the i-th direction of Ag + A and 
df Ag + A for the upper boundary. In particular, x belongs to d^Ag + A if it belongs 
to Ai + a and ±xt = I + A. Since the measure ha m ,k is uniform, 

B w M , K [(^9)s'] = E^ m iK [{T v g)g'] 

if the support of T x g and the one of r y g do not intersect the one of g' . Therefore, 
choosing A large enough, the previous sum vanishes. This proves that the function 
(|6.5(l has mean zero with respect to all canonical invariant measures. 

At this point, we follow the classical approach of nongradient systems (cf. 0, 
Chapter 7) to estimate the second term of (|6.4|) using the standard Rayleigh- 
Schroedinger perturbation theorem for the largest eigenvalue of a symmetric oper- 
ator. After a few steps we bound the difference of the second term of (|6.4(l with 
eN 2 ~ d D N (f) by 



N 2 ~ d e x - 
M d ^ 



In this formula, B stands for the function l|6.5l) . the supremum is carried over all 
integers < K < \Am\ and < •, • >^ A K is the inner product in L 2 (/ja m _k)- Since 
the spectral gap of the generator of the symmetric exclusion process in Am is of 
order M 2 and M 2 N~ 1 vanishes as TV f oo, by the perturbation theorem for the 
largest eigenvalue of a symmetric operator, the previous expression is less than or 
equal to 

1 ^suv<{-L% m )- 1 B,B>^ m>k . 
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Consider the linear term in the variational formula for the H-\ norm of B. It is 
given by 2 < B,f >/j Aju k f° r some function / in L (/m m ,k)- Since B has mean 
zero with respect to all canonical invariant measures, this is in fact a covariance 
that we estimate by Cq(£)M 2 + C\M~ 2 < /, / > PA K . By the spectral gap for the 
symmetric exclusion process, the second term is bounded by < (—L A )f, f > MA K 
if we choose C sufficiently small. Therefore, < (— L s Au )~ x B, B > MA K is bounded 
by C(£)M 2 . Since we are in dimension d > 3, the last displayed equation van- 
ishes as TV j co. This proves that the second term in (|6.4|) may be estimated by 
eN 2 ^ d Dj^(f) and an expression which vanishes as N f oo. □ 

We have just proved that we may replace L* by L A in (j4.3|l . We show now that 
we can replace the average ?y M (0) by the average i] e (0). 

Lemma 6.2. Fix a function f in 5, a smooth function G : R_|_ x T d — * R and 
assume that M satisfies the conditions 13.2(1 . For every T > 0, 

pT 



lim limsup 

£-►00 jv->oo 



dt / N 1 - 4 J2 G{t,z/N) 



, { L* At (A e f) (n M (0) , r,) - L* Af (A £ f) [rf (0) , r,) } f t 



= 



Proof. We have seen in the proof of the previous theorem that it is enough to show 
that 

N^ d J2 G^z/N^Ll^f^iO)^)-^^)^^),^} 
is negligible. 

Consider a class of function B(/3,rj), < j3 < 1, whose support is contained 
in Af. Repeating the well known steps of the proof of the one block estimate we 
obtain that 

|A M I . 

B( V M (0),r,)f( V )dv» = £Cfc(/) B(K/\A M \,t,)f MtK (rj)dn* 

K=0 

where, 

Jm 



IA M ,K 



C kU) = / 1{ 2^ r,{x)=K}fdv» , f MiK (rj) = j 

and /m is the conditional expectation E v n [f | Tm ] ■ Here, for a set A, T A stands 
for the tr-algebra generated by {r](z), z £ A}. At this point, B(K/\Am\, ■) is a local 
function with support in Ai and we repeat the procedure for fu,K, MAm.-K" m pl ace 
of /, v£ . We obtain in this way that the previous sum is equal to 

I A M | |A<| 

^2 C K (f)^2C k (f M ,K) / B(K/\A M \,r))f M ,K,i,k(v)dVA t ,k 

K=0 fe=0 J 

with the obvious definitions for Ck{fM,n), fM,K,e,k- 

Using that the Dirichlet form is convex, we may estimate 

f N 1 ~ d £ G(t,z/N)B(T, M (p),r,)(T- g f)dv5 - eN 2 - d D N (f) 
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by 

I A M | |A f 

N ~ d E E C*(/*)£ WXr,*) (6-6) 

zST ^ K=0 k=0 

{G(t,z/N)N J B(K/\A M \, r,)f Z M,K A Md» Al , k - ^D Ae (f^ KAk ^ Ae , k )} . 

In this formula, f z — T— Z f and D A( (•, fiA £ .k) is the Dirichlet form associated to the 
generator L S A( and the reversible measure HK t ,h- Assume that B(K/\Am\, rf) has 
mean zero with respect to all invariant states fi\ e ^, which is the case of the func- 
tion we are considering in this lemma. By the Rayleigh-Schroedinger perturbation 
theorem for the largest eigenvalue of a symmetric operator, the expression inside 
braces in the previous formula is less than or equal to 

< (-Lir'm/^MlvlBiK/l&Mlr)) >„ e , k ■ (6.7) 

We claim that in the particular case of this lemma, the previous expression is 
bounded by Ce~ l {K l\A M \-k/\Ai\) 2 . Indeed, let h be the local function f(K/\A M \, 
rf) — f(k/\Ag\, rf). In the case where B is the function which appears in the statement 
of the lemma, the linear term of the variational formula for the H—± norm is 

rr~T E / ( L * T v h )f d ^e,k > 

where / is in L 2 (/jA^,fc)- Since L*r y h is a local function which has mean zero with 
respect to all invariant measures, we may localize / around y, replace the scalar 
product by a covariance, use the spectral gap of the symmetric exclusion process, 
restricted to a cube whose length depend only on the support of h, and apply 
Schwarz inequality to bound < (VbE[f\ T/v]) 2 > by < (Vb/) 2 >. At the end we 
obtain that the previous expression is less than or equal to 

( ^ < (L*Tyh) 2 > MAt . + <-L S f,f>„.. . 



Since f(-,?y) is smooth, uniformly in rj, L*r y h is absolutely bounded by |AT/|Am| — 
fc/|A|| |. This proves that is bounded above by Ce- 1 (K/\A M \ ~k/\A e \) 2 . 

Up to this point we proved that the expression inside braces in (jfi.fijl is bounded 
above by Cs~ 1 (K/\Am\ — k/\Ai\) 2 . Recalling the definition of the constants ap- 
pearing in (|6.6|l . we have that this sum is in fact 

C 



eN d 



{r? M (0)-^(0)} (r_,/)di£. 



It remains to apply the two blocks estimate to conclude the proof. 



6.3. Replacement of Hij (77) by <TijF(r] e (0)) + Ji,j(v (ty)- Fix a smooth function 
G : T d x M. + — > R and two function f, g in Since the local functions f(/3, •) have 
a common finite support, for each fixed y, there exists a finite integer A such that 

z£A A 
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Since f are smooth functions, the difference between the previous expression 
and 

z£A A 

is absolutely bounded by C(A,f) |?y M (0) — ^(0)|, for some finite constant C(A,f). 
By the two blocks estimate, the average over Tjy of this absolute value is negligible. 
After this replacement, the third line of i|4.3[) is seen to be composed of three 
different types of terms: 

£ G(t,x/N)r x £ P *(y)y i y j r,(m-V(y)} , 

x£f d N y£l d 

£ G(t, x/N) r x £ p*(y) Ul „(Q)[1 - „(!,)] r#(„) , 

x£l d N y£l d 

where, for some function f) in 5", 

- V 0) , £ f)(^(0),r z ry). 

By the one block estimate, the first sum can be replaced by 
£ G{t,xlN)o i4 F{rf{x)) . 

x£T% 

We claim that the second sum is negligible because 77(0) [1 — v(y)]^f nas mean 
zero with respect to all canonical invariant measures. Indeed, repeating the steps 
of the one block estimate, we are reduced to estimate 

sup / r)(0)[l-r)(y)]V o , v £ f(K/\Ai\,T z r)) dfi Ae ,K , 

z£A A 

where the supremum is carried over all < K < |A^|. A change of variables 
£ = (T°' v 7] permits to rewrite the previous expression as 

sup / {n{y) - ri(0)} £ f(K/\Ae\,T z r))dnA e , K ■ 

z£Aa 

The integral vanishes for each fixed K because [iA e ,K is a uniform measure. 

The third type of term requires some notation. For a function ()(/?, 77), smooth 
in the first coordinate and with a common finite support in the second, let 

t)(a,(3) = E Vf3 [\){a,rj)] . 

For 1 < i, j < d and y in Z d , let 

ttfifrv) = v(P)[i-v(v)]Vo,v £ U(P,r.v)Vo,v £ fMw)- 

z£A A zeA A 

Notice that f) is smooth in the first coordinate and have a common finite support 
on the second coordinate. Moreover, an elementary computation shows that 
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In this formula, < ■, ■ >p stands for the inner product in L (yp). Denote the right 
hand side by Jf^fJp). Lemma ffi . 31 below shows that we may replace in ()4.3|) the 
third type of term by 

]T G{t,x/N)J Mi {rf{x)). 
Up to this point, we proved that the third line of H4.3[) is equal to 

d 

(!/ 2 ) E E (d u Mt,x/N){d Uj \){t,x/N){o^Ftf(x)) + J n ,ijtf(x))} 

plus a term of order o(N d ), where 

Jn^tf) = 2J2 <kn(/3,-),(-i s )r,f,,„(/3,-)>/3 ■ 

xGZ d 

Recall the definition of the function Ji y j(f3) given in 14. 6[) . By Theorem 15.11 with 
the notation introduced in section [21 the previous sum is equal to 

d 

(V2) E E (duMt^/NKdujXK^x/N^Firfix)) + Jij(tf{x))} 
».i=i xer d N 

plus a term of order o(f, N d ). To conclude this subsection, it remains to prove the 
next result. 

Lemma 6.3. Fix a function \)(P,rj) smooth in the first coordinate and with finite 
common support in the second. For positive integers I, m, let 

Vljrj) = -L J2 f)(^(0),r^)-t)(^(0),^(0)) 



2/6A„ 



Then. 



lim hmsupsupf f -L V T^Jjf) f dv" - eN 2 - d D N (f)\ = 
/or all e > 0. 

Proof. Since v-y is the Bernoulli product measure, for each (5, an elementary com- 
putation shows that (c> 7 )f)(/3, 7) = XLeA < h(f3, rj);r}{x) > 7 , where < ■; • > 7 stands 
for the covariance with respect to v n and A for a finite set which contains the 
common support of the function f)(/3, •). In particular, the derivative (e? 7 f))(/3, 7) is 
uniformly bounded. Hence, 

t,(^(0),^(0))-^(r/(0),rr(0))| < C(&) |r? m (0) - rf{0) | 

for some finite constant C(f)). It follows from the two blocks estimate that we may 
replace f)(i/(D), i/(0)) by 6(^(0), ?7 m (0)) in the definition of V% m . 

Following the classical proof of the one block estimate, we are reduced to estimate 

SU P / TT1 E f)(^/|Ad,^)-f)(^/|A,U m (0)) dn Ai , K , 

where the supremum is carried over all < K < |A^|. For each fixed £, denote 
by Kt the integer which maximizes the previous variational formula. There exists 
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a subsequence I' such that Ke> /\Ae>\ converges to sonic density j3 in [0,1]. In 
particular, the limsup, as i ] oo, of the previous expression is less than or equal to 

yeA m 



sup 

i9e[Q,i] 



|Ar, 



because the finite marginals of the canonical measure converges to the grand canon- 
ical measures. Since fj(/3, •) is a smooth function, 

In particular, the previous variational formula is bounded above by 
sup / TlTl bi^Tyf])- E^[t)(/3, V )} dvp + C sup / \ri m (0)-p\dvp . 

This expression vanishes as m f oo because vp is a product measure and f)(/3, •) 
are local functions with a finite common support. This concludes the proof of the 
lemma. □ 



6.4. Estimation of the current. Fix i < i < d and recall the definition of V i (rj) 
given just after JO}. Let 



A Mf (t, V ) = N 1 -* Y (^ui A) (t, x/N)t x \f (r)) 



By the nongradient estimates, for every T > 0, 

limsup limsup / dt v% (d-q) (rj) A itN/A (t,r)) 



< Co sup 

a£[0,l] 



i<j<d 

for some finite constant Cq. Here ||| • ||| is the norm introduced at the beginning of 
section [5J We refer to section 6 of jS] for the proof. Note that we don't need in 
the present context the multiscale analysis of [H]- By Theorem 15 . 1 1 this expression 
vanishes if we replace fi by fi !rt and let n f 00. 
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